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We present the general relativistic calculation of the energy release associated with a 
first order phase transition (PT) at the center of a rotating neutron star (NS). The en- 
ergy release, -Ercli is equal to the difference in mass-energies between the initial (normal) 
phase configuration and the final configuration containing a superdense matter core, 
assuming constant total baryon number and the angular momentum. The calculations 
^ ' are performed with the use of precise pseudo-spectral 2-D numerical code; the poly- 

QQ ' tropic equations of state (EOS) as well as realistic EOSs (Skyrme interactions. Mean 

, Field Theory kaon condensate) are used. The results are obtained for a broad range of 

■ metastability of initial configuration and size of the new superdense phase core in the fi- 

nal configuration. For a fixed "overpressure" , 5P, defined as the relative excess of central 
pressure of a collapsing metastable star over the pressure of the equilibrium first-order 
PT, the energy release up to numerical accuracy does not depend on the stellar angular 
momentum and coincides with that for nonrotating stars with the same 5P. When the 
equatorial radius of the superdense phase core is much smaller than the equatorial ra- 
dius of the star, analytical expressions for the iS^ei can be obtained: Ej-ei is proportional 
to (SP)'^-^ for small SP. At higher SP, the results of 1-D calculations of Ej.gi(5P) for 
non-rotating stars reproduce with very high precision exact 2-D results for fast-rotating 
' stars. The energy release-angular momentum independence for a given overpressure holds 

$H ' also for the so-called "strong" PTs (that destabilise the star against the axi-symmetric 

perturbations), as well as for PTs with "jumping" over the energy barrier. 

Keywords: dense matter — equation of state — stars: neutron - stars: rotation 



1. Introduction 

Many theories of dense matter predict that at some density larger than the nu- 
clear saturation density, a phase transition (PT) to some "exotic" state (i.e. boson 
condensate or quark deconfinement) occurs; for review see e.g., [1-3]). A first-order 
PTs are particularly interesting from the astrophysical and observational point of 
view, because are associated with a meta-stable state of dense matter. One can 
thus expect, in the case of PTs occurring in the interior of NSs, the release of 
non-negligible amount of energy. Here we will focus on the basic features of the 
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Fig. 1. Left panel: Strong and weak PTs on mass-radius diagram (dotted line marks unstable 
configurations). Right panel; the energy release S^el "vs the over-pressure parameter hP for strong 
(two upper curves) and weak (lowest curve) PTs. Points are colored differently for different total 
angular momenta J = (0, 0.1, 1.3) X GMq/c and follow the curve for the J = (non-rotating) 
configurations. 



energy release-angular momentum independence; for complete description we refer 
the reader to [4, 5]. The text is arranged as follows: in Sect. [2] we briefly describe 
the results of calculations. Sect. [3] contains conclusions and open questions. 



2. Calculation of the energy release 

The hydrostatic, axi-symmetric and rigidly rotating configurations of compact stars 
with and without PTs were obtained using the numerical GR library LORENE 
( http : //www ■ lorene . obspm . f r] ) . We assume that the baryon mass and the 
total angular momentum J are constant during the PT from the metastable con- 
figuration C to stable configuration C* . From the microscopic point of view the 
first-order PT is characterised by the over-pressure parameter 5P = [Pc — Po)/Pq, 
where Pc is the central pressure of the configuration C and Pq is the equilibrium pres- 
sure, at which the PT occurs. We distinguish between the so-called weak and strong 
PTs, characterised by the density-jump parameter A — p* /p, where p and p* are the 
densities of normal and condensed phase at Pq. PTs with A > ^{1 + Pq/ pc^) (strong 
PTs) destabilise the stars with arbitrarily small cores against the axi-symmetric 
perturbations [6-8] . The energy release (energy difference) is defined as 

E,,i = c^[M{C)^MiC*)]M^j (1) 

and it is quite remarkable that (up to numerical accuracy) it does not depend on 
the rotation state of the configuration i.e. on the total angular momentum J of the 
star and agrees well with the E^ei for non-rotating stars. In Fig. [1] we present the 
results for the polytropes (parametric EOSs in the form of P = Kri^ , where n is the 
number baryon density, K is is the pressure coefficient and F is called the adiabatic 
index; detailed description of parameters used can be found in [5] ) . 
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Fig. 2. Left panel: The strong PT with the "jumping" over the energy barrier on the mass-radius 
diagram. Right; Energy release E^gi vs the over-pressure parameter SP for the strong PT. Points 
with different colors correspond to PTs with different angular momentum J = (0, 0.1, . . . , 0.9) X 
GMq/c^ also for negative 5P. Solid line - result for J = (non-rotating stars). Point A corresponds 



to Co 



■ Cg , point B to Co 



■ Cq' = Co, point C to C2 > CJ and point D to Cn 



In the case of realistic EOSs the results are qualitatively the same. We have 
checked the phenomenon for different realistic EOSs i.e. Skyrme interactions or the 
kaon condensate, presented as an example here. The crust obeys the EOS of Douchin 
& Haensel [9]. The matter below the PT is described using the relativistic mean- 
field theory, [10]. The dense phase is the kaon condensate, with coupling of kaons 
to nucleons proposed by Glendcnning & Schaffncr-Biclich [11], the optical potential 
C/^" being equal to —115 MeV. The stability of configurations with Pc below that 
for the equilibrium PT is particularly interesting. If the star is excited initially, 
e.g. is pulsating, then the formation of a large dense phase core is possible, but it 
requires climbing ("jumping") over the energy barrier associated with formation of 
a small core. The results are shown in Fig. [2j 



3. Conclusions and open questions 

Our numerical calculations show the E^eiiJ) independence during the weak as well 
as strong first-order PTs, for large rotation rates and large oblatnesses of the stars: 
it is therefore not an property of slow rotating stars only. The independence holds 
also for PTs with negative over-pressure SP i.e. when the configuration "jumps" 
over the energy barrier to reach another stable configuration. For small positive dP 
analytical relations were found: E^rci oc (SP)^'^. The energy release Ei-d ^ 10^^ — 10^^ 
erg is an absolute upper bound on the energy which can released in such PT. In 
astrophysical situation, the energy can be distributed between stellar pulsations, 
gravitational radiation, heating of stellar interior. X-ray emission from the neutron 
star surface, and even a gamma-ray burst. 

Currently there is no mathematical proof of the J-indcpcndcncc of E^-d. It may 
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be interesting to look at the problem from the "thermodynamical" point of view. 
We write the total energy of the star (gravitational mass M) as 



2nj + 2 [ PNy^cfix, 



(2) 



where ^ is the baryon chemical potential, u* the time component of the fluid 4- 
velocity, A the number of baryons in the star, the angular velocity, P the fluid 
pressure, N the lapse function and y^cfix the covariant volume element in the 
constant t hypersurface Sj. Eq. ([2]) is valid for any axi-symmetric stationary and 
rigidly rotating fluid star which obeys a barotropic EOS, as established by Bardeen 
& Wagoner in 1971 [12]. Since C and C* have the same baryon number A and the 
same angular momentum J, we get the energy release 
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The terms on the right-hand-side of Eq. [3] are of comparable magnitude in their 
contribution to the energy release, so the possible "smallness" of some of them 
in relation to the others is not responsible for the E,-e\{J) independence. We will 
address the problem of mathematically proving this property in the near future. 
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